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Introductory Remarks 


The application of mathematical nonlinear programming algorithms has enjoyed 
considerable success in the automated structural synthesis environment. A 
sizeable majority of research pertaining to optimum structural design has 
primarily focussed on statically loaded structures. The addition of dynamic 
loading in the form of harmonic excitation (Ref. 1) or random loads (Ref. 2) to 
the structure introduces significant new problems. Since all aerospace vehicles 
operate in a dynamic loads environment, it is imperative that efficient 
automated structural sizing procedures be developed that adequately account for 
these complexities. Furthermore, an increasing focus on multidisciplinary 
design and optimization dictates a realistic representation of each 
participating discipline. 

The optimum sizing of airframe structures is performed in a complex analysis 
environment where structural stability must be carefully monitored in addition 
to elastic deformation under static and dynamic loads. Hence an analysis tool 
that accounts for the interaction of structural deformations and unsteady 
airloads, is a necessary ingredient to the development of any synthesis 
capability. The present report documents the framework of an optimization 
system, developed around the ISAC computer program (Ref. 3), with SPAR (Ref. 4) 
and CONMIN (Ref. 5) as the structural analysis and optimization programs, 
respectively. The NOS control language feature on the CDC Cyber computer is 
used to direct the flow between various segments of the system. In its present 
form, the system retains a degree of modularity, permitting substitutions of the 
analysis packages. 

The thrust of the present work was directed towards establishing an optimization 
capability for sizing airframe structures that are subjected to a combination of 
deterministic and random flight loads. The random vibration environment 
introduces the need for selecting a statistical process that best describes the 
random loads and permits computation of the dynamic response parameters of 
interest. Furthermore, it requires a formulation of design constraints that 
would minimize the conservativeness in the design and retain computational 
viability. In the present work, the random loads are treated as a stationary, 
homogeneous process with a Gaussian probability distribution. A frequency 
domain analysis is selected for the solution of the dynamic response problem. 
This permits a power spectral density input for the gust velocities, a 
representation that is typically available from experimental observations. The 
ISAC system of programs has the feature to incorporate aribitrary gust velocity 
input spectra. 

Subsequent sections of this report discuss the formulation of the analysis 
problem, the structure of the optimization programming system and a 
representative numerical example. An appendix with procedure files and FORTRAN 
pre- and post-processor listings is included for completeness. 







Dynamic Response Analysis 


In a frequency domain analysis, the system equations of motion can be written as 

( -JW + CK] - [A]) {W} = {G} (1) 

where [M] and [K] are the mass and stiffness matrices, respectively; {W} is the 
displacement vector; [A] is the matrix of loads due to oscillation at frequency 
(0, and {G} is the force coefficient array per unit gust velocity. As is typical 
of large structural configurations, a modal approach is adopted wherein a finite 
number of elastic modes are used to model the structural deformations. The 
displacement {W} can be represented approximately by the superposition of the 
characteristic modes 


{W} = [<D] {5} (2) 

where 1^1 is the matrix of 'm' normalized eigenmodes and {5I is an m-dimensional 
vector of generalized coordinates. Substituting equation (2) into equation (1) 
and then transforming the resulting equation from the frequency domain into the 
Laplace domain yields the equations of motion in the following form 

[M^^s^ + (1 + - A^^]{5} = {5} (3) 

where M is the generalized mass matrix; A and G* are the generalized 

aerodynamic force matrices due to motion and gust respectively; g is the 
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structural damping coefficient and is the natural frequency associated with 

the 5-th generalized coordinate. Equation 3 is a system of linear, 
nonhomogenous equations that are solved for the vector of generalized 
coordinates. A similar approach can be adopted to include a control surface, 
and the option to do so exists in the ISAC system. 

Typical response parameters for which design constraints can be formulated 
include stresses, nodal displacements, and accelerations at specified 
locations. The focus of the present work was in implementing structural 
reliability constraints for which element stresses are the primary response 
quantities of interest. In the modal displacement approach, the stress at a 
location j is computed as 


where c.. is the stress coefficient at location 'j' due to a unit displacement 

' J 

in the i-th generalized coordinate 5^. For a random gust input, these quanti- 
ties are obtained as root mean square (rms) values which are then used in the 
constraint formulation. 


Structural Reliability Constraints 

For a structure subjected to random gust loads, failure can be caused by a 
single overstress or as a result of cumulative damage in fatigue. These failure 
conditions are termed by Johnson (Ref. 6) as single excursion and fatigue. 
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respectively. Constraints must be formulated for both these failures modes. 

For a stress response S(t) that is stationary and has a Gaussian probability 
, distribution, an assumption that large values of S(t) arrive independently, 

leads to a_Poisson probability function for the number of times *n' that a large 
magnitude S is exceeded in time 't'. If the occurrence of this large stress 
♦ magnitude is not permitted in time T , the single excursion constraint, g<-^r» 

can be written as ^ 








exp ( *- ) - 1 < 0 

n C 


(5) 
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Fatigue damage can be estimated on the basis of the Palmgren-Miner theory and 
results in a fatigue constraint, of the following form 


9c = 



( 2 ) 


b-2 b-1 

-y^s 


r 1 < 0 


( 6 ) 


where, is the specified fatigue life; r is the gamma function; 0 ^ and 

are the rms values of the stress and stress rates and b and c are constants 
obtained from empirical relations for different materials. 

Failure constraints derived from a fracture mechanics standpoint have also been 
proposed (Ref. 7) where the constraints expressed in equations (5) and (6) are 
essentially unified into a single constraint. The initiation and propagation of 
cracks in a material results in both a reduced allowable stress and can also be 
viewed as accumulated damage in the structure. In the present work constraint 
equations (5) and (6) were implemented in the design. 


Optimization System 

This section of this report gives a brief description of the optimization system 
and its component programs. The optimization system combines a finite element 
program SPAR, a general-purpose, nonlinear programming- based optimization 
program COMMIN, and a system of aeroelastic response analysis programs ISAC. 

The flow between these programs is controlled in a manner similar to the PROSSS 
system for optimum design (Ref. 8). The role of each of these programs is 
discussed next. The reader is referred to Appendix A for a listing of the 
various routines described in this section. 

SPAR 


The program is a system of independent processors which communicate with each 
other through a data base which consists of libraries of named data sets. The 
structural model geometry definition, material specification and element connec- 
tivity are defined in a SPAR input runstream NREPT, and the data set library is 
stored as NRLA. This is the non- repeatable part of the structural analysis pro- 
gram and does not have to be re-executed as the program iterates in the optimi- 
zation mode. The SPAR program is used to generate for the ISAC system, quanti- 
ties such as mode-shapes, eigenvalues, generalized masses and stress coeffi- 
cients defined in equation (4). It is also used to generate finite-difference 
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gradients for the eigenmodes and frequencies. A FORTRAN program GRAD generates 
an input runstream corresponding to the set of new design variables. This input 
runstream, SECPROP, is executed before the repeatable segment REPT is called to 
obtain the new response parameters. A FORTRAN program EIGSTR employs the SPAR 
data handling utilities (Ref. 9) to extract the desired information from the 
data set libraries and store the information for post processing. 


I SAC 


The I SAC system of routines is employed to compute the aeroelastic response of 
the flight vehicle. The structural node deflections corresponding to the vibra- 
tion modes obtained in SPAR are related to the deflections and slopes on aerody- 
namic boxes by a 2-D, spline interpolation technique in DLIN. These aerodynamic 
slopes and deflections are used by a doublet lattice program DLAT, in the compu- 
tation of the generalized aerodynamic forces for a range of reduced 
frequencies. The generalized mass, stiffness and aerodynamic matrices are 
passed to the routine DYNARES where a gust spectrum is selected and the aero- 
elastic response parameters are obtained. For a random gust input, the rms 
responses can be computed by numerical integration of the response spectra 
within this routine. In the present organization, however, DYNARES writes a 
state space representation of the system equations and passes the system 
matrices to another group of programs PADLOCS (Ref. 10), where the rms responses 
are computed. Additionally, this program provides correlation coefficients 
between various response components which are required in the formulation of 
stress constraints in the event of a combined stress failure criterion 
(Ref. 11). The input files required for ISAC execution are listed in 
Appendix A. For a detailed description of the input parameters, the reader is 
referred to Reference 3. 

CONMIN 


CONMIN is a general purpose optimization subroutine for solving constrained and 
unconstrained, linear or nonlinear, continuous problems. It is based on the 
method of feasible-usable search directions (Ref. 5) and was adapted in the 
present framework in a piecewise linear mode of operation. The general optimi- 
zation problem for CONMIN can be posed as follows 
Minimize a function 

F(d^, d2...dj^) (7) 

Subject to the inequality constraints 

gj(d^, d 2 ,...d^) 40 j = 1,2, ...m (8) 

and the side constraints on the design variables d^. 

d/s< d. 4 dj“ i = l,2,...n (9) 

In general, both the objective function and constraints can be nonlinearly 
dependent on the design variables. In a piecewise linear mode of operation, the 
objective function and design constraints are regarded to be linear functions of 
the design variables, eliminating the need to compute function derivatives at 
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every iteration for a general nonlinear problem. A linear approximation to a 
general nonlinear problem requires that constraints be imposed on design vari- 
able change for a given cycle of optimization. An appropriate choice of design 
variables (reciprocal variables) will often enhance the linearity in the design 
space. 


Program Organization 

A flow chart depicting the implementation of the optimization system is shown in 
Figure 1. The corresponding procedure file PRLOOP is listed in Appendix A. The 
structural model is initialized for the starting design and the eigenmodes, 
natural frequencies and response coefficients and their sensitivity to perturba- 
tions in design variables are computed. This information is then used to com- 
pute the response constraints sensitivity in a finite difference operation. The 
two sensitivity computations were separated on purpose. If the eigenvalue and 
eigenvector sensitivities are assumed constant, the outer gradient loop can be 
executed once in the optimization sequence. The inner gradient loop for con- 
straint sensitivities is the only required operation for subsequent cycles of 
optimization. By the same token, if the eigenvalue/eigenmode sensitivities are 
known to be more nearly constant over a wider variation in design variables, the 
outer loop can be executed a lesser number of times in the optimization 
sequence. The objective function, constraints and their gradients are passed to 
an optimization program where the design variables are altered within the limits 
imposed by the side constraints. This sequence is repeated until convergence. 


Numerical Examples 

A built up finite element model of a wing structure was developed as a test 
problem for the program implementation described in the preceeding sections. 

The cross-section and planform of this stringer-membrane, cantilever structure 
is shown in Figures 2 and 3. The model element sizes were preprocessed by a 
fully-stresssed design cycle for 1-g cruise loads. These member sizes were 
established as the lower bounds for the gust design. A Dryden spectrum with a 
gust intensity of 2.54 m/sec (typical of storm conditions) was specified in the 
input. The first six elastic modes were chosen to represent the structural 
deformation and were considered adequate for such a cantilever structure. In 
more realistic designs for gust loads, rigid body modes in plunge and pitch 
should also be included. The addition of these in the present system consti- 
tutes no added complexity. The design variables selected for this minimum 
weight design were the cross-sectional areas of the bar elements. An allowable 
stress of 1761.4 kg/cm^ was used in the computation of the first excursion 
failure constraints. First excursion and fatigue failure lives of 2 years were 
specified in the design process. The design converged in 13 piecewise linear 
cycles to the values shown in Table 1. The first cantilever mode in bending 
dominates the stress distribution in the structure, as is evidenced by a concen- 
tration of material at the root section. The first excursion constraint is 
active at the optimum. 


Concluding Remarks 

This report documents the development and implementation of an optimization 
system around ISAC for the structural design of airframe structures with aero- 
elastic constraints. The programming system is regarded as the necessary first 
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step for further studies in the systematic design of aircraft structures, 
operating in a random gust load environment. Although developed primarily for 
the design of structures under random gust loads, the program in its present 
form is suitable for inclusion of static stress, displacement and frequency 
constraints. An s-plane formulation of the flutter constraint can also be added 
with minor modifications. These features coupled to a built-in control analysis 
capability, make this programming system an extremely powerful tool for multi- 
disciplinary analysis and design. Furthermore, the programming system provides 
a natural test-bed for multilevel optimization studies. The computational 
efficiency of the system can be significantly enhanced with the inclusion of an 
analytical gradient capability. 
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Figure 1. Flowchart depicting the organization of the optimization system 
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membrane elements 



Figure 2. Cross section of wing structural box 




Figure 3, 


Finite element model of Wing Structural Box. Node numbers nl , n2 represent 
lower and upper surface co-ordinates, respectively. Panel elements were 
sized by a fully-stressed design approach. 


































Design 

Variable 

Element Connectivity 

Initial Design 
(cm^) 

Minimum Gage 
(cm^) 

Final Design 
(cm2) 

1 

1-9,9-17,8-16,16-24 

7-15,15-23,6-14,14-22 

10.3226 

10.1096 

17.9516 

2 

17-25,25-33,24-32,32-40 

22-30,30-38,23-31,31-39 

10.3226 

4.6297 

9.4458 

3 

33.41,41.49,40-48,48-56 

39-47,47-55,38-46,46-54 

10.3226 

0.6452 

2.5239 

4 

2-10, 10-18, 3-n ,11-19 
4-12,12-20,5-13,13-21 

10.3226 

10.1096 

17.9574 

5 

18,26,26-34,19-27,27-35 

20-28,28-36,21-29,29-37 

10.3226 

4.6309 

9.5703 

6 

34-42,42-50,35-43,43-51 

36-44,44-52,37-45,45-53 

4.516 

0.6452 

2.5298 

7 

54- 62,62-70,70-78 

55- 63,63-71 ,71-79 

56- 64,64-72,72-80 
49-57,57-65,65-73 

4.516 

0.6452 

1 .1761 

8 

53-61 ,61-69,69-77 
52-60,60-68,68-76 
51-59,59-67,67-75 
50-58,58-66,66-74 

4.516 

0.6452 

0.9661 

Objective 

Function Weight 

286.33 Kg 

195.56 

270.41 Kg 


Table 1 


Numerical results for the cantilevered wing structural box 








Appendix A 


LIST,FN=SUBFIL 

PRLaOF,TibOOO. 

USER , 358625C , PATHA. 

CHARGE , 1 02632 , LRC . 

DELIVER. 1229 PH PRQC TEST CASE 
* 

* THIS FILE IS PRLQOP 

* 

MAP, OFF. 

BEGIN.PRLOOP. 

CGPYSBF , I NPUT , OUTPUT . 

DA YF I LE , L=DAYPRLP . 

REPLACE , DAYPRLP . 

EXIT. 

CDPYSBF , I NPUT , OUTPUT . 

D A YF I LE , L=D AYPRLP . 

REPLACE, DAYPRLP. 

EQI ENCOUNTERED. 

/ 


Above procedure file is used to initiate a batch run for optimization 



PROCFIL is a listing of all procedure files used in the present 
optimization system. 


LIST,FN=PROCFIL 

-PRDC,DYNA. 

GET,DCM=DCMV2A/UN=BS7010N. 

BET , D YNARES=DYNAV3A/UN=aa70 1 ON . 

REVERT. 

.PROCjCONT. 

GET , CONM I NB/UN=753437N . 

GET,CONB. 

COP YBR , CONB , PRAT . 

COP YBF , CONM I NB , PRAT . 

GET , DEB I GN , OLDOB J , CONMDT , CONP AR . 
REWIND,*. 

PRAT. 

REWIND,TAPE6. 

APPEND , CONOUT , TAPE6 . 

REPLACE, OLDOB J. 

REPLACE, DESIGN. 

REVERT. 

.PROC,OBFUN. 

RETURN,*. 

GET , OB J FUB , DEC I S , I TER , DEB I GN , C VD 1 , OBD J . 
PACK,CVD1. 

OBJFUB. 

REPLACE, CONMDT. 

RETURN,*. 

REVERT. 

.PROC,DYNDES. 

6ET,TAP2G. 

BEGIN, DYNA, 

BEGIN, PRDES. 

BEGIN, DYNA. 

BEG IN, PRDES. 

BEGIN, DYNA. 

BEGIN, PRDES. 

BEGIN, DYNA. 

BEG IN, PRDES. 

BEG IN, DYNA. 

BEG IN, PRDES. 

BEGIN, DYNA. 

BEG IN, PRDES. 

BEG IN, DYNA. 

BEG IN, PRDES. 

BEG IN, DYNA. 

BEG IN, PRDES. 

BEGIN, DYNA. 

BEG IN, PRDES. , 

REVERT. 

. PROC,PRGUST. 

GET,TAPE5=P0UT. 



BEGIN,DCM. 

DCM, , , , ,DCMIN1. 

RETURN , DCM , TAPES . 

REWIND,TAPE9. 

REPLACE , T APE9=PHT9 . 

BEGIN, DLIN. 

GET ,TAPE1=W2P0S1 , TAPE3=TAPE31 , TAPES. 
DLIN. 

REWIND,*. 

RETURN , TAPE 1 , TAPES . 

GET , TAPE 1 =W2P0S2 , TAPE3=TAPE32 . 

DLIN. 

REWIND,*. 

RETURN , TAPE 1 , TAPES . 

GET , T APE 1 =W2P0S3 , T APE3=T APE33 . 

DLIN. 

REWIND,*. 

REPLACE , TAPE9=PHT9 . 

RETURN,*. 

REVERT. 

.DATA,DCMIN1. 

BATCH 

STORE 

GMASS -1 1 324 1 
F 

GENERALIZED MASSES 
STORE 

FRED -1 1 18 1 
F 

FREQUENCIES IN HZ 
STORE 

LOADS -1 1 1296 1 
F 

STRESS COEFFICIENTS 
STORE 

DAMPINGS -1 1 18 1 
F 

DAMPING COEFFICIENTS 
STORE 

GMASS -1 2 324 1 
F 

GENERALIZED MASSES 
STORE 

FREQ —1 2 18 1 
F 

FREQUENCIES IN HZ 
STORE 

LOADS -1 2 1296 1 
F 

STRESS COEFFICIENTS 
STORE 

DAMPINGS -1 2 IB 1 
F 

DAMPING CONSTANTS 
STORE 
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□MASS -1 3 324 1 
F 

GENERALIZED MASSES 
STORE 

FREQ -1 3 18 1 
F 

FREQUENCIES IN HZ 
STORE 

LOADS -1 3 1296 1 
F 

STRESS COEFFICIENTS 
STORE 

DAMPINGS -1 3 18 1 
F 

DAMPING CONSTANTS 
END 

.PROC,PRDES. 

GET,TAPE9=PHT9. 

COPYBR , TAP2G , TAPE2 . 

REWIND,TAPE2. 

DYNARES. 

REPLACE , TAPE 1 =T APE7 . 

RETURN, *,TAP2G. 

GET , DESDAT , TAPE? , DES I GN=DEZ IGN. 

DESIGN, DESDAT. 

REWIND,CVDATA. 

GET, REDUB. 

REDUB. 

APPEND, CVDl ,CVDATA. 

RETURN, »,TAP2G. 

REVERT. 

.PRDC,PREXT. 

GET , EXTPLB,DECIS, DESIGN, CONDAT, ITER, NODAL. 
REWIND,*. 

EXTPLB. 

REWIND,*. 

REPLACE , WP0S=W2P0S 1 , WPOS 1 =W2P0S2 , WP0S2=W2P0S3 . 
REPLACE, POUT. 

REVERT. 

.PROC, STORE. 

GET , E I GSTRB , I TER , QRD AT A , LOO 1 =RRL A . 

REWIND,*. 

El GSTRB. 

REPLACE, GRDATA. 

REVERT. 

.PRQC,PRLOaP. . 

SET(R2=10) 

SET(R3=0) 

lOPT.CTIME. 

SET<R1=0) 

MAP , OFF . 

WHILE, R1.LT.R2,F INI. 

SET(R1=R1+1) 

DISPLAY (Rl) 

GET , GRADE , I TER , DEC I S , GRDATA , OBD J . 

GRADE. 
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REPLACE , I TER , CONDAT , DBD J . 
RETURN, GRDATA. 
IFE,R1.LT.R2,JNP1. 
BEBIN,PRRSPAR. 

BEGIN, STORE. 

RETURN,*. 

ENDW,FIN1. 

ENDIF, JMPl. 

BEGIN, PREXT. 

BEGIN,PRGUST. 

GET , T APE5=T APE5T , T APE9=PHT9 . 
BEGIN,DLAT. 

ALPHA. 

REPLACE , T APE9=PHT9 . 

RETURN,*. 

BEGIN,DYNDES. 

BEGIN,OBFUN. 

GET,CVD1=CV. 

REPLACE, CVDl. 

SET(R3=R3+1) 

BEGIN, CONT. 

GET , DEC I S , DEB I GN , I TER , PQSB . 
POSB. 

REPLACE, DECIS. 
IFE,R3.LT.9,JMP2. 
IFE,FILE(IFLG,.NOT.AS) , JMP3. 
RETURN,*. 

GOTO, lOPT. 

ENDIF, JMP2. 

ENDIF, J MP3. 

REVERT. < — 

. PRDC , PRCSTR , X X X=W2ZP0S . 

* 


PRLOOP is the procedure file 
to execute the various segments 
of the optimization sequence. 


* EXECUTE THIS PROCEDURE AS FOLLOWS: 

* -DLNPROC,DLNPROC. 

* 

RETURN, LGO. 

GET,DLNLIB. 

*PUT HERE A GET TO YOUR TAPES - INPUT FILE 

*NEXT MAKE AVAILABLE A SPAR FILE WHICH TO EXTRACT DATA FROM 
GET,L001=RRLA. 

BET,X=EIGSTR. 

SET, GRDATA, ITER. 

REW I ND , X , I TER , GRDATA . 

FTN , I =X , L=0 , 0PT=0 , ST AT I C . 

GTR , LGO , LL . REL/E I GSTR 
GTR , DLNL I B , LL . REL/BLKD AT . 

RENAME, LGO=LL. 

LDSET ( L I B=DLNL I B , MAP=SBEX /OUT ) 

LOAD , L.u]0 . 

NDGO(EIBSTRB) . 

REVERT . 

.PROC,PRRSPAR. 

. * 

. * REPEATABLE SPAR 



GET ,SPAR=SPAR14I ,DCU=DCU14I/UN=750756N 
GET , SPARLA=NRLA. 

REWIND, SECPROP. 

SPAR , SECPROP , SPOUT . 

RETURN, SPOUT. 

REWIND,*. 

GET,REPT. 

SPAR , REPT , SPOUT . 

REWIND, SPARE A. 

REPLACE , SP ARLA=RRLA . 

RETURN,*. 

REVERT. 

-PROC,PRNSPAR. 

. * 

.* NONREPEATABLE SPAR 

. * 

GET , SPAR=SPAR 141, DCU=DCU 141 /UN=750756N 
GET,NREPT. 

SPAR , NREPT , SPOUT . 

REPLACE , SP ARLA=NRL A . 

RETURN,*. 

REVERT. 

.PROC,DLIN. 

GET , DCM=DCMV2A/UN=aa7010N. 

GET , DL I N=DL I NV3A/UN=8G70 1 ON . 

REVERT. 

.PROC,DLAT. 

GET , DCM=DCMV2A/UN=aa70 1 ON . 

GET , ALPHA=DLATV5A/UN=aa70 1 ON . 

REVERT. 

.PROC,DCM. 

GET , DCM=DCtiV2A/UN=8870 1 ON. 

REVERT. 

.PROC,PLOT. 

REWIND.SAVPLT. 

PLOT . TEKPOST , TE ( AUTO ( O ) ) . 

NOTE. /:G507> 

REVERT. WHAT A CRUMMY PLOT! ! 
.PROC,PLTRVAX. 

REWIND, S A VPLT. 

PLOT. TEKPOST, TE (AUTO (O) .RVAX) . 
N0TE./:G5O7> 

REVERT. 

EOI ENCOUNTERED. 

/ 
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LIST,FN=CONT 

PROGRAM CONOPT ( CONPAR , DLDOB J , DBS I GN , CONMDT , I FLG , 

♦T APE56=CDNP AR , TAPE 1 2=DLD0B J , 

*T APE3 1 =DES I GN , T APE22=C0NMDT , TAPE 1 4= I FLG , OUTPUT , 

*TAPE6) 

C 

C *♦**■»(•*■»(•***•«■***#**•«■#**♦♦********************#*****•*♦****■*!•♦****** 

C 

C THIS IS THE OPTIMIZER WITH PIECEWISE LINEAR APPROXIMATION 
C 

C PRABHAT HAJELA JUNE 1985 
C 

C *♦**##■»!■•)(•******#***#******•«•**♦****•)<■***********■«■****»*******#**** 

C 

COMMON/CNMNl /DELFUN , DABFUN , FDCH , FDCHM , CT , CTMIN , CTL , CTLMIN , 
♦ALPHAX , ABOB J 1 , THETA , OBJ , NDV , NCON , NS I DE ,T PR I NT , NFDG , NSCAL , L I NOB J 
*ITMAX, ITRM,ICNDIR,IGOTO,NAC, INFO, INFOG, ITER 
C0MM0N/CNMN2/RDUM (50) , IDUM (50) 

C0MM0N/C0NSAV/REAL(50) ,INT(25) 

DIMENSION X (50) , VLB (50) ,VUB(50) ,G(200) ,SCAL(50) , 

*DF(50) , A (50,200) ,S(50) ,G1 (200) ,G2(200) ,B (200,200) , 

#C (200) , ISC (200) , IC (200) ,MS1 (400) 

DIMENSION VLBTMP(50) ,VUBTMP(50) ,XINIT(50) ,GI (200) , 

*GRDG (50,200) jGRDOBJ (50) 

DATA NDV/8/,NC0N/144/ 


COMPUTE CONMIN DIMENSIONS 

NLIM1=NDV+1 

Nl=NDV+2 

N2=NC0N+2*NDV 

N3=N2+5 

N4=MAX0(N3,NDV) 

N5=2*N4 


READ CONMIN PARAMETERS 

READ ( 56 , * ) N I TER , I PR I NT , NS I DE , NFDG , L I NOB J , I TMAX , I CND I R , NSCAL , 

* I TRM , I GOTO , NDV , NCON , FDCH , FDCHM , CT , CTM I N , CTL , CTLM I N , THETA , PH I , 
*DELFUN , DABFUN , ALPHAX , ABOB J 1 

READ OLD OBJECTIVE FUNCTION VALUES 

READ(12,*)0BJ1 ,0BJ2,0BJ3,T0L 
C 

DO 5 1=1, NCON 
5 ISC(I)=0 

C 

C PRESCRIBE BOUNDS ON DESIGN VARIABLES FOR OVERALL OPTIMIZATION 
C 
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VLBTMP(1)=1.567 
VLBTMP (2) =0.7176 
VLBTMP(3)=0. 1 
VLBTMP(4)=l-567 
VLBTMP (5) =0.7178 
VLBTMP(6)=0. 1 
VLBTMP (7) =0. 1 
VLBTMP (8) =0. 1 
DO 6 1=1 ,NDV 
VUBTMP(I)=5.0 
6 CONTINUE 

C 

C READ INITIAL DESIGN VARIABLES 
C 

READ(31 ,*) <XINIT<I) ,1=1 ,NDV) 

PRESCRIBE BOUNDS FOR LOCAL OPTIMIZATION CYCLE 

DO 12 1=1, NDV 
VUB(I)=1.4*XINIT<I) 

VLB(I)=0.6*XINIT(I) 

X (I)=XINIT(I) 

IF ( VUB ( I ) . GT. VUBTMP ( I ) ) VUB ( I ) =VUBTMP ( I ) 

IF (VLB < I) .LT. VLBTMP (I> ) VLB ( I ) =VLBTMP < I > 
CONTINUE 

READ INITIAL FUNCTION AND GRADIENT INFORMATION 

READ(22,#) (GI (I) ,I=1,NC0N) ,0BJ1 
*, ( (GRDGd , J) , 1=1 , NDV) , J=1 ,NCON) , 

♦(GRDOBJ (I) ,1=1 ,NDV) 

OBJ=OBJl 
lO CONTINUE 

CALL TO CONMIN OPTIMIZER 

CALL CONMIN (X, VLB, VUB, G,SCAL,DF, A, S, Gl ,G2,B,C, 
*ISC, IC,MS1 ,N1 ,N2,N3,N4,N5) 

IF(IG0TD.ED.0)60 TO 555 
IF ( INFO. EQ. 2) GO TO 24 

FUNCTION AND CONSTRAINT VALUES 

SUM=OB J 1 
DO 23 1=1, NDV 
DX=X (I)-XINIT(I) 

SUM=SUM+GRDOBJ ( I ) *DX 
; CONTINUE 

OBJ=SUM 

DO 22 J=1,NC0N 
SUM=GI (J) 

DO 21 1=1, NDV 

DX=X(I) - XINIT(I) 

21 SUM=SUM+GRDG(I,J)*DX*(XINIT(I)/X(I) ) 

G(J) =SUM 

22 CONTINUE 
GO TO 100 

24 CONTINUE 
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c 

C GRADIENT INFORMATION 
C 

DO 25 1=1, NDV 

25 DF ( I ) =GRDOBJ ( I ) 

NAC=0 

DO 26 J=l,NCON 

IF(G(J) -LT.CT)GO TO 26 

NAC=NAC+1 

IC(NAC)=J 

26 CONTINUE 

DO 27 1=1, NDV 
DO 28 J=1 ,NAC 
J1=IC<J) 

A(I ,J) =GRDG(I ,J1) 

28 CONTINUE 

27 CONTINUE 
C 

GO TO 100 
C 

C TERMINATE OPTIMIZATION LOOP 
C 

555 CONTINUE 

REWIND 31 

WRITE (31 ,*) (X (I) , 1=1 ,NDV) 

C 

C TEST FOR CONVERGENCE OVER THREE CYCLES 
C 

0BJ3=0BJ2 
0BJ2=0BJ1 
0BJ1=0BJ 
REWIND 12 

WR ITE ( 1 2 , * ) OBJ 1 , OB J2 , OB J3 , TOL 
IF (ABS (0BJ2) -LE. 0.0001) GO TO 666 
DA=ABS ( (DBJ3-0BJ2) /0BJ2) 

DB=ABS( (0BJ2-0BJ1) /0BJ2) 

C 

C TERMINATE PROGRAM 
C 

IF (DA. GT. TOL. OR. DB. GT. TOL) GOTO 666 
WRITE(14,*)NC0N 
666 CONTINUE 
STOP 
END 

EQI ENCOUNTERED. 

/ 
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LIST,FN=eiGSTR 

PROGRAM E I GSTR ( GRD ATA , I TER , TAPE2= I TER , T APE32=GRD ATA , 
♦OUTPUT) 

C 

C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

c 

C THIS PROGRAM READS AND STORES THE DATA REQUIRED FOR 
C THE EIGENANALYSIS AND STRESS COEFFICIENT GRADIENTS 
C 

C PRABHAT HAJELA JUNE 1985 
C 

C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

c 

DIMENSION G1 (100) ,G2(13000) ,G3(10000) ,64(10) ,X (10) ,Y(10) 
DIMENSION AM0DE(3,B0) ,FREQ(25) ,STRS(72,6, 1) ,STR1 (10000) 
DATA NL I B/ 1 / , NMODE/6/ , ND4 1/10/, NSTR/72/ , P I /3 . 1 4 1 59/ 

♦ , NNODE/SO/ , NELT/72/ , NDOF/240/ 

C 

C READ FREQUENCY 
C 

CALL DAL(NLIB, 11 ,FREQ,0,0,KADL,IERR,NW,NE,LB, IT, 

♦4H V I BR , 4HEVAL , 4HMASK , 4HMASK ) 

C 

C STRESS COEFFICIENTS 
C 

CALL DAL (NLIB,11,STR1(1) ,0,0,KADL,IERR,NW,NE,LB,IJ, 

1 4HSTR3,4HE23 ,1,1) 

CALL DAL (NLIB,11 ,STR1 (433) ,0,0,KADL, IERR,NW,NE,LB,IT, 

1 4HSTRS,4HE23 ,1,2) 

CALL DAL (NLIB,11,STR1 (865) ,0,0,KADL,IERR,NW,NE,LB,IT, 

1 4HSTRS,4HE23 ,1,3) 

CALL DAL (NLIB,11,STR1 (1297) ,0,0,KADL,IERR,NW,NE,LB,IT, 

1 4HSTRS,4HE23 ,1,4) 

CALL DAL (NLIB,11,STR1 (1729) ,0,0,KADL,IERR,NW,NE,LB,IT, 

1 4HSTRS,4HE23 ,1,5) 

CALL DAL (NLIB,11,STR1 (2161) ,0,0,KADL,IERR,NW,NE,LB,IT, 

1 4HSTRS,4HE23 ,1,6) 

C 

C OBJECTIVE 
C 

CALL DAL(NLIB, 11 ,0BJP,0,0,KADL,IERR,NW,NE,LB, IT, 

♦4H0B JF , 4HAUS , 4HMASK , 4HMASK ) 

C 

C EIGENMODES 
C 

CALL DAL (NLIB,10,AMQDE,0,0,KADR,IERR,NWDS,NE,LB,IT, 

1 4HV I BR , 4HM0DE , 4HMASK , 4HMASK) 

C 

C 

READ ( 2 , ♦ ) I TER , ND V 

NCDEF=NELT^NMQDE 

NLIM1=NDV+1 

NLI M2=NMGDE*NL IM 1 

NL I M3=NC0EF*NL I M 1 

NL I M 4=ND0F^NL I M 1 ♦NMODE 
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c 

C IF ITER.EQ.O CREATE GRDATA 
C 

IF(ITER.GT.O)GO TO 2000 
DO 1000 1=1,NLIM2 

1000 G1(I)=0.0 

DO 1001 I=1,NLIM4 

1001 G2(I)=0.0 

DO 1002 I=1,NLIM3 

1002 G3(I)=0.0 

DO 1003 I=1,NLIM1 
G4(I)=0.0 

1003 CONTINUE 

WRITE (32,*) (G1 (I) ,I=1,NLIM2) , 

*(G2(I) , 1=1 ,NLIM4) , 

*(G3(I) ,I=1,NLIM3) , 

*(G4(I) , 1=1 ,NLIM1) 

2000 CONTINUE 
C 

REWIND 32 

C IF ITER.GT.O READ DATA FROM OLD GRDATA FILE 
C 

READ(32,*> (G1 (I) ,I=1,NLIM2) , 

*<62 (I) , 1=1 ,NLIM4) , 

*<63(I) ,I=1,NLIM3) , 

*<G4(I) ,I=1,NLIM1) 

REWIND 32 
C 

C STORE DATA 
C 

C OBJECTIVE 
C 

G4(ITER+1)=0BJP 

C 

C FREQUENCY 
C 

DO 2 I=1,NM0DE 
FREQCY=SQRT (FREQ ( I ) ) / (2. *PI ) 

K= ( ITER*NMODE) +I 
61 (K)=FREQCY 
2 CONTINUE 

C 

C EIGENMODES 
C 

I D= I TER*NDOF»NMODE 

DO 4 I=1,NM0DE 

KS=K ADR + ( I - 1 ) *NSECTS ( LB ) 

CALL RI0(NLIB,2,2,KS,AM0DE,LB) 

DO 6 J=1,NN0DE 
DO 7 K=l,3 
ID=ID+1 

G2(ID)=AMQDE(K, J) 

7 CONTINUE 

6 CONTINUE 

4 CONTINUE 

C 

C STRESS COEFFICIENTS 
C 
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ID=ITER*NCOEF 

JC=0 

DO 8 I = l,h4MODE 

IDl=a-l)*NELT»6 

DO 9 J=1,NELT 

ID2=(J-1)*6 + 5 + IDl 

IDK=ID2+1 

ID=ID+1 

G3(ID)=STR1 (IDK) 

9 CONTINUE 

CONTINUE 

WRITE OBJECTIVE, El BEN- INFO AND STRESS COEFFICIENTS 

WRITE (32,*) (G1 (I) , 1=1 ,NLIN2) , 

*(G2(I) ,I=1,NLIM4) , 

*(G3(I) ,1=1 ,NLIM3) , 

*(G4(I) , 1=1 ,NLIM1) 

REWIND 32 
STOP 
END 

EOI ENCOUNTERED. 

/ 
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LIST,FN=EXTPLT 

PROGRAM EXTPLT ( DEC I G , DEB I GN , CDNDAT , POUT , NODAL , I TER , 
*WPQS , WPOS 1 , WPDS2 , 

*TAPE33=DES I GN , TAPE22=C0NDAT , TAPE7=P0UT , TAPEl =NODAL , 
♦TAPE4=WP0S1 ,TAPE8=WP0S2, 

*T APE3 1 =DEC I S , T APE9= I TER , T APE3= WPOS , I NPUT , OUTPUT ) 

C 

C *****************************»■*■***********•»(■********#*•»••* 

C 

C THIS PROGRAM PROVIDES A NEW SET OF EXTRAPOLATED MODES 
C EIGENFREQUENCIES AND RESPONSE COEFFICIENTS TO ANALYSIS 
C 

C PRABHAT HAJELA JUNE 1985 
C 

C ********************************************************* 

C 

DIMENSION B1 (100) ,G2<1500) ,G3(1100) ,G4(10) ,DD(60) 
DIMENSION GRDG(8, 1500) ,GRDG1 (8,1100) ,GRDG2(8, 10) , 
♦GRDOBJdO) 

DIMENSION NODE (80) ,X (10) ,Y(10) ,XX (10) ,AMODE (6,3,80) 
DIMENSION G ( 1000) , FREQ (10) ,DELT( 10) ,VECN( 1500) 
DIMENSION GMAS(54,54) ,S1 (100) ,S2(10000) ,STRS (72,54, 1 ) 
DATA NMODE/6/ ,ND0F/240/ ,NELT/72/ ,NE41/10/ , 

*NSN0DE/40/ ,NN0DE/80/ 

OPEN ( UN I T=7 , F I LE= ' POUT ' ) 

OPEN ( UN I T= 1 , F I LE= ' NOD AL ' ) 

OPEN (UNIT=33,FILE= ' DESIGN ' ) 

OPEN ( UN I T=22 , F I LE= ' CONDAT ' ) 

OPEN ( UN I T=3 1 , F I LE= ' DEC I S ' ) 

OPEN ( UN I T=9 , F I LE= ' I TER ' ) 

OPEN ( UN I T=3 , F I LE= ' WPOS ' ) 

OPEN (UNIT=4,FILE=' WPOS 1 ' ) 

OPEN ( UN I T=8 , F I LE= ' WP032 ' ) 

READ(9,*) ITER,NDV 
NCOEF=NELT *NM0DE 
N1=NDDF*NM0DE 
NLIM1=NDV+1 
NL I M2=NM0DE*NL I M 1 
NL I M3=NC0EF*NL I M 1 
NLIM4=ND0F*NLIM1 
NSTR=NELT 
DO 40 1=1,18 
DO 41 j=i,ia 
GMASd , J)=0.0 
IFd.EQ. J)GMAS(I,J)=1.0 
41 CONTINUE 

40 CONTINUE 

READ(31 ,*) (X (I) , 1 = 1 ,NDV) , (Yd) , 1 = 1 ,NE41) ,FDCH 
READd,*) (NODE (I) ,I = 1,NSN0DE) 

REhD(22,») (G1 (I) , 1=1 ,NMODE) , (G2 ( I ) , 1=1 ,N1 ) , 

*(G3(I) ,I=1,NC0EF) ,G4(1) , 

»(GRD03J (I) , 1=1 ,NDV) , 

(GRDGd ,J) ,1 = 1 ,NDV) ,J=1 ,N1) , 

* ( (GRDG2( I , J) , 1=1 ,NDV) , J=1 ,NMODE) , 

* ( (GRDGl (I , J) , 1=1 ,NDV) , J=1 .NCOEF) 
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DO 30 INN=1,NLIM1 

READ (33,*) (XX (I) ,I=1,NDV) 

REWIND 33 

. _IF(INN. EQ. IJ GQ JP.31 _ _ 

XX (INN-i)=)CX(iNN-i) *('iI+FDCH) 

1 CONTINUE 

UPDATE STRESS COEFFICIENTS 

DO 2 I=1,NC0EF 
STR=G3 ( I ) 

DO 4 J=1,NDV 
DELT(J)=XX(J)-X(J) 

STR=STR + GRDGl (J,I)*DELT(J) 
CONTINUE 
G ( I ) =STR 
CONTINUE 

FREQUENCIES 

DO 6 I=1,NM0DE 
EIG=G1(I) 

DO 8 J=1,NDV 
DELT(J)=XX(J)-X<J) 

E I G=E I G+GRDG2 ( J , I ) *DELT < J ) 

CONTINUE 

FREQ(I)=EI6 

CONTINUE 

EIGENMODES 

DO 10 1=1, N1 
VEC=G2 ( I ) 

DO 12 J=1,NDV 
DELT(J)=XX(J)-X(J) 

VEC= VEC+GRDG ( J , I ) *DELT ( J ) 

12 CONTINUE 

VECN ( I ) =VEC 
10 CONTINUE 

ID=0 

DO 14 I=1,NM0DE 
DO 16 J=1,NN0DE 
DO 17 K=l,3 
ID=ID+1 

ANODE ( I , K , J ) = VECN (ID) 

17 CONTINUE 

16 CONTINUE 

14 CONTINUE 

DO 21 K=1,NM0DE 
DO 20 I=1,NSN0DE 
IF(INN.LE.3) THEN 

WRITE (3, 1000) (ANODE (K,J,NODE(I) ) ,J 
1 OOO FORNAX ( 3E 1 5 . a ) 
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GO TO 19 
EMDIF 

IF ( INN. LE. 6) THEN 

WRITE<4, 1000) (AMODEIK, J,NODE(I) ) ,J=1 ,3) 
GO TO 19 
END IF 

IF ( INN. LE. 9) THEN 

WRITE<8, 1000) (AMODEIK, J,NODE(I) ) ,J=1 ,3) 
ENDIF 

19 CONTINUE 

20 CONTINUE 

21 CONTINUE 

IJ= ( INN-1 ) *NMODE 
DO 50 IK=1,NM0DE 
IJK=IJ+IK 

50 SKI JK ) =FREQ ( I K ) 

IL=(INN-1) *NCOEF 
DO 51 IM=1,NC0EF 
ILM=IL+IM 

51 S2(ILM)=G(IM) 

30 CONTINUE 

JP=0 

DO 52 J=1,NLIM2 
DO 53 I=1,NELT 
DO 54 K=l,l 
JP=JP+1 

STRS(I,J,K)=S2(JP) 

54 CONTINUE 

53 CONTINUE 

52 CONTINUE 
IC=NLIM2/18 
DO 61 IZ=1,54 

61 DD(IZ)=0.01 

DO 60 1=1, IC 

WRITE (7,*) ( (GMASdl , J) ,J=1 , IS) , 1 1 = 1 , 18) 
KP=(I-1)*18 

WRITE(7,*) (SI (KP+II) ,11=1,18) 

DO 55 IEL=1,NELT 
DO 56 K=1 , 1 

WRITE (7,^^^) (STRS ( lEL, IY,K) , IY=KP+1 ,KP+18) 
56 CONTINUE 

55 CONTINUE 

WRITE (7,*) (DD(IZ) ,IZ=1 ,18) 

60 CONTINUE 

STOP 
END 

EOI ENCOUNTERED. 

/ 
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LIST,FN=GRAD 

PROGRAM GRAD ( DEC I S , I TER , GRD AT A , COND AT , SECFROP , FL AG , DBD J , 
*TAPE3 1 =DEC I S , TAPE30= I TER , T APE32=GRDATA , TAPE22=C0NDAT , 

»T APE7=SECPROP , TAPE 1 4=FL AG , T APE2=0BD J ) 

C 

C *******#****#************»***#***************»#****#*** 

c 

C THIS PROGRAM IS A PREPROCESSOR TO A GRADIENT COMPUTING 
C SEGMENT - GRADIENTS OF EIGENVALUES, EIGENVECTORS AND 
C STRESS COEFFICIENTS ARE COMPUTED 
C 

C PRABHAT HAJELA, JUNE 1985 
C 

C *******************************#**»-tHHHH(-********-IHf-***** 

C 

DIMENSION GRDG<B, 1500) ,GRDDBJ (10) ,X (10) ,G1 (100) 
*,62(13000) ,G3(10000) ,G4(10) ,Y(10) 

*,6RDG1 (8,1100) ,GRDG2(8,10) 

NM0DE=6 
ND41=10 
NSTR=72 
ND0F=240 
NCOEF=NSTR*NMODE 
0PEN(UNIT=30,FILE=' ITER' ) 

OPEN (UNIT=31 ,FILE= ' DECIS ' ) 

OPEN ( UN I T=7 , F I LE= ' SECPROP ' ) 

OPEN ( UN I T=22 , F I LE= ' COND AT ' ) 

RE AD ( 30 , * ) I TER , ND V 
REWIND 30 

READ DESIGN VARIABLES AND FINITE DIFFERENCE STEP SIZE 

READ (31,*) (X(I) , I = 1,NDV) , (Yd) ,I = 1,ND41) ,FDCH 
ITER=ITER+1 


IF ITER.EQ.NDV+1 THEN ALL DV'S PERTURBED 

IF(ITER.EQ.NDV+1)G0 TO lOO 
IFdTER.GE. 1) X (ITER)=X(ITER)*(1.+FDCH) 

WRITE SPAR INPUT RUNSTREAM 

REWIND 7 
WRITE(7,9S) 

98 FORMAT (' CXQT TAB',/, 'BO 
DO 1 1=1, NDV 
WRITE(7,*) I ,X (I) 

1 CONTINUE 

WRITE (7,99) 

99 FORMAT ('SA') 

DO 2 1=1 ,ND41 
WRITE(7,*) I ,Y(I) 
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2 CONTINUE 

WRITE(7, 1010) 

1010 FORMAK'CXQT EXIT') 

REWIND 7 

WR I TE ( 30 , * ) I TER , NDV 

GOTO no 

100 CONTINUE 

COMPUTE GRADIENTS 
ITER=-1 

WRITE (30,*) I TER, NDV 
NLIM1=NDV+1 
NL I M2=NM0DE*NL I M 1 
NL I M3=NC0EF*NL I M 1 
NL I M4=ND0F*NL I M 1 *NMQDE 
OPEN (UNI T=32 , F I LE= ' GRD AT A ' ) 

READ (32,*) (G1 ( I ) , 1=1 ,NLIM2) , 

*(G2(I) ,1=1 ,NLIM4) , 

*(G3(I) , 1=1 ,NLIM3) , 

*(G4(I) , 1=1 ,NLIM1) 

OBJECTIVE 

DO 10 1=1, NDV 
J1=I+1 

GRDOBJ ( I ) = (G4 ( J 1 ) -G4 ( 1 ) ) / (FDCH*X ( I ) ) 
CONTINUE 

FREQUENCY 

DO 25 J=1,NDV 
II=NMODE*J 
DO 20 I=1,NM0DE 
II1=II+I 

GRDG2 (J,I) = (G1(II1)-G1(I))/ (FDCH*X ( J) ) 
CONTINUE 
CONTINUE 

STRESS COEFFICIENTS 

DO 27 J=1 ,NDV 
II=MCOEF*J 
DO 22 1=1 ,NCOEF 
II1=II+I 

GRDGl (J,I) = (G3(II1) -G3 ( I ) ) / (FDCH*X ( J ) ) 
CONTINUE 
CONTINUE 

EIGENMODES 

N1=ND0F*NM0DE 
DO 29 J=1,NDV 
I I=N1*J 
DO 24 1=1, N1 
II1=II+I 

GRDG ( J , I ) = (G2 (I I 1 ) -G2 ( I ) ) / (FDCH*X ( J) ) 
CONTINUE 
CONTINUE 
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WRITE GRAD INFD 

WRITE (22,*) (G1 (I) ,I=1,NM0DE) , 

*(G2(I) ,I=1,N1) , 

*(G3(I) , 1=1 ,NCQEF) ,G4(1) , 

*(GRDOBJ ( I ) , 1=1 ,NDV) , 

*( (GRDGd ,J) ,I = 1,NDV) ,J=1,N1) , 

*((GRD62(I,J) ,I=1,NDV) ,J=1,NM0DE) , 

*( (GRD61 (I , J) , 1=1 ,NDV) , J=1 ,NCOEF) 

REWIND 22 

OPEN ( UN I T=2 , F I LE= ' 03D J ' ) 

WRITE(2,*)G4(1) , (GRDOBJ (I) , 1=1 ,NDV) 

EXIT FROM GRADIENT LOOP BY DECLARING TAPE 14 AS LOCAL 

0PEN(UNIT=14,FILE='FLAG' ) 

WRITE(14,*) ITER 
110 CONTINUE 
STOP 
END 

EOI ENCOUNTERED. 

/ 
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LIST,FN=OBJFUN 

PROGRAM OB JFUN ( DEC I S , I TER , DES I GN , CVD 1 , CONMDT , OBD J , 
*TAPE8=CVD 1 , TAPE 1 0=DES I GN , TAPE 1 6=C0NMDT , 

*TAPE4=DEC I S , TAPE7= I TER , TAPE2=0BD J ) 

C 

c 

C THIS PROGRAM COMPUTES AND WRITES THE SENSITIVITY 
C INFORMATION AS REQUIRED BY THE OPTIMIZATION 
C PROGRAM, CONOPT. 

C 

C PRABHAT HAJELA JULY 1985 
C 

CC *-»f-»****#****-»***##*#*****#*****-»****************** 

C 

DIMENSION A (12000) ,STR(1000) ,STRR(1000) 

DIMENSION Y(IO) ,X(10) ,XX<10) ,TEM(1000> 

DIMENSION GSE(IOOO) ,GFF ( 1000) ,GRDOBJ ( 10) ,GRDG ( 10 ,200) 
DATA NSTR/72/,SIGAL/25000. /,TSE/6. 1516SE+7/ , PI /3. 14159/, 
*TFF/6. 15168E+7/ ,GAM5/24. / ,CEE/1 . E+41/ ,ND41/10/ 
OPEN(UNIT=a,FILE='CVDl ' ) 

OPEN (UNIT=10,FILE= 'DESIGN ' ) 

OPEN (UNIT=16,FILE= 'CONMDT ' ) 

QPEN(UNIT=4,FILE='DECIS' ) 

OPEN <UNIT=7,FILE=' ITER' ) 

OPEN ( UN I T=2 , F I LE= ' OBD J ' ) 

READ (7,*) ITER,NDV 

READ(2,*)DBJ1 , (GRDOBJ(I) ,I=1,NDV) 

READ (10,*) (X (I) , 1=1 ,NDV) 

READ (4,*) (XX (I) ,I=1,NDV) , (Y(I) ,I=1,ND41) ,FDCH 

NLIM1=NDV+1 

NL I M2=NL I M 1 *NSTR 

NSTR1=1 + 2*NSTR 

NSTRSO=NSTR 1 **2 

NDIM= (NSTRSQ+NSTRl ) /2 

READ(8,*)N1 

DO 2 IJ=1,NLIM1 

READ(8,*)N1, (Ad) ,I = 1,NDIM) 

IN=(IJ-1)*NSTR 

JCK=1 

DO 4 I=1,NSTR1-1 
TEM(I)=SQRT(A(JCK) ) 

JCK=JCK+NSTR1- ( I-l ) 

4 CONTINUE 

DO 6 I=1,NSTR 
K=NSTR+I 

STR(I+IN)=TEM(I) 

6 3TRR(I+IN)=TEM(K) 

2 CONTINUE 

C 

C CONSTRAINT - SINGLE EXCURSION FAILURE 
C 
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DD 10 I=1,NLIM2 
STR<I)=STR(I)*100. 

STRR ( I ) =STRR ( I ) *100. 

FUN= ( S I GAL**2) / ( 2 . *STR ( I ) *STR ( I ) ) 

FUN=-FUN 

IF (FUN. LE. -675.0) FUN=-675. O 
I F ( FUN . GE . 74 1 . 0 ) FUN=74 1 . 0 

10 GSE(I)=(TSE/PI)*(STRR(I)/STR(I> )*EXP(FUN) - 1. 

C 

C CONSTRAINT - FATIGUE FAILURE 
C 

DD 12 I=1,NLIM2 
FUN=(TFF/ (PI*CEE) )*8.*6AM5 
12 GFF(I)=(FUN*STRR(I)*(STR(I)**7) ) - 1. 

C 

C CONSTRAINT DERIVATIVES 
C 

DO 14 1=1, NDV 

K=I*NSTR 

DO 16 J=1,NSTR 

GRDG ( I , J ) = (GSE (K+J) -GSE ( J ) ) / (FDCH*X ( I ) ) 

16 CONTINUE 

14 CONTINUE 

DO 18 1=1, NDV 
K=I*NSTR 

DO 20 J=NSTR+1 ,2*NSTR 

GRDG ( I , J) = (GFF (K+J-NSTR) -GFF ( J-NSTR) ) / (FDCH*X ( I ) ) 
20 CONTINUE 

15 CONTINUE 
C 

C WRITE OBJECTIVE, CONSTRAINTS AND THEIR DERIVATIVES 
C 

WRITE (16,*) (GSE (I) , 1=1 ,NSTR) , 

* (GFF ( I ) , 1=1 ,NSTR) ,0BJ1 , 

*( (GRDG(I,J) ,1=1, NDV) ,J=1,2*NSTR) , 

*(6RD0BJ(I) ,1=1, NDV) 

C 

STOP 

END 

EOI ENCOUNTERED. 

/ 
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LIST,FN=POST 

PROGRAM POST ( DES I GN , DEC I S , I TER , T APE2= I TER , 

*T APE7=DEC I S , TAPE8=DES I GN ) 

C 

C *****#*•«•****************************■)«•*•»*******♦•»■♦*•*(■******** 

c 

C THIS PROGRAM UPDATES THE DEC IS FILE AND CAN ALSO BE USED 
C TO DETERMINE IF MODES NEED TO BE RECOMPUTED IN PROBLEM 
C 

C PRABHAT HA J EL A JULY 1985 
C 

C *•***#***♦************#■•*•******•»***•«■**•«•*■«•*•*•******■•(■*■»********* 

c 

DIMENSION X(10) ,Y(10) ,XX<10) 

DATA ND41/10/ 

OPEN (UNIT=2,FILE=' ITER • ) 

OPEN ( UN I T=7 , F I LE= ' DEC IS') 

OPEN ( UN I T=8 , F I LE= ' DES I GN ' ) 

READ(2,*) ITER,NDV 

READ(7,*) (XX(I) ,I=1,NDV> , (Y<I) ,I=1,ND41) ,FDCH 
READ(8,*) (X(I) ,I=1,NDV) 

DO 2 1=1, NDV 
2 XX(I)=X<I) 

REWIND 8 
REWIND 7 

WRITE<7,100) <XX(I) ,1=1, NDV) , (Yd) ,I=1,ND41) ,FDCH 
100 FORMAT (FIO. 4) 

STOP 

END 

EOI ENCOUNTERED. 

/ 


^ ’ ™RolRftM REDUC < CVDflTA , TAPE3-CVDATA , OUTPUT > 

c ««.»...*....»...♦.»**********•»»*♦******•***•*************** 

C THIS PROGRAM REDUCES THE SIZE OF THE COVARIANCE MATRIX 

c 

cc 

c PRABHAT HAJELA JULY 1985 

C ************************************************************ 

C 

DIMENSION CV(145,145) 

OPEN (UNIT=3,FILE= 'CVDATA ' ) 

READ(3,+i^)Nl 

READ ( 3 ,^) < ( CV (I,J) ,J=1,NI) ,1-1, Nl) 

REWIND 3 
WRITE(3,»)N1 
DO 2 1=1, Nl 

WR I TE < 3 , * ) ( CV ( I , J ) , J= I , N 1) 

2 CONTINUE 

CL0SE(UNIT=3) 

STOP 

END 

EDI ENCOUNTERED. 


A-21 



LIST, 

FN=NREPT 


CXQT 

TAB 


START 80,4 5 6 


MATS 

1 10 

. 5-*“6 .3 .1 


2 10 

. 5‘+*6 .3 .1 


JLOC 

1 

0. 00 

-67.50 


0.00 

-67.50 

"T 

0.00 

-37.50 

4 

0.00 

-7.50 

5 

0. 00 

22.50 

6 

0. 00 

22.50 

7 

0.00 

-7.50 

a 

0.00 

-37.50 

9 

40. 00 

-81.57 

10 

40 . 00 

-81.57 

1 1 

40 . 00 

-53.61 

12 

40.00 

-25.64 

13 

40. 00 

2.32 

14 

40. 00 

jC— m 

15 

40.00 

-25.64 

16 

40 . 00 

-53.61 

17 

80.00 

-95.64 

18 

80.00 

-95.64 

19 

80 . 00 

-69.71 

20 

80.00 

-43. 79 

21 

80 . 00 

-17.86 

22 

80 . 00 

-17.86 

23 

80. 00 

-43. 79 

24 

SO. 00 

-69.71 

25 

1 20 . OO 

-109.71 

26 

1 20 . 00 

-109.71 

27 

120.00 

-85.82 

28 

120. 00 

-61.93 

29 

120. 00 

-38.04 

30 

120.00 

-38.04 

31 

120. 00 

-61.93 

“TO 

120.00 

-85.82 

• 

1 60. 00 

-123. 78 

34 

160. 00 

-123.78 

35 

160. 00 

-101.92 

3<£j 

160. 00 

-80.07 

37 

1 60 . 00 

-58.22 

33 

1 60 . 00 

-58.22 

39 

1 60 . 00 

-80. 07 

40 

160. 00 

-101.92 

41 

200.00 

-137.84 

42 

200.00 

-137. S4 

43 

200.00 

-1 18.03 

44 

200 . 00 

-93. 22 

45 

200. 00 

-78. 40 

46 

200 . 00 

-78. 40 


non-repeatable part of SPAR 


3.78 
7.02 

10.80 

10.80 
7.02 
3.78 
0.00 
0. 00 
3.89 
6.91 
10.43 
10.43 
6.91 
3.89 
0.37 
0.37 
4.00 

6.80 
10.07 
10.07 

6.80 
4. 00 
0.73 

0. 73 
4.11 
6. 69 
9. 70 
9.70 
6.69 
4.11 

1 . 10 
1 . 10 
4.22 
6.58 

9.33 
9 . -Jisil' 
6. 58 
4 . 22 

1.47 

1 . 47 
4. 33 

6. 47 
B.97 

а. 97 

б. 47 

4.33 
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47 

200 . 00 

-98.22 

1.83 


48 

200 . 00 

-113.03 

1.83 


49 

240.00 

-151.91 

4.44 


50 

240.00 

-151.91 

6 . 36 


51 

240.00 

-134. 14 

8 . 60 


52 

240.00 

-116.36 

8.60 


53 

240. UO 

- cm 

• wJO 

6. 3o 


54 

240.00 

-93.53 

4.44 


55 

240. 00 

-116.36 

2.20 


56 

240.00 

-134. 14 

2. 20 


57 

280.00 

-165.98 

4.55 


58 

230.00 

-165.98 

6.25 


59 

280 . 00 

-150.24 

8.23 


60 

280 . 00 

-134.50 

8.23 


61 

230 . OO 

-118. 76 

6.25 


62 

280.00 

-118.76 

4.55 


63 

280.00 

-134.50 

2.57 


64 

230.00 

-150.24 

2.57 


65 

320.00 

-180.05 

4.66 


66 

320.00 

-180.05 

6. 14 


67 

320. 00 

-166.35 

7.87 


68 

320. 00 

-152.64 

7.87 


69 

320.00 

-138. 94 

6. 14 


70 

320. OO 

-138.94 

4.66 


71 

320.00 

-152.64 

2.93 


72 

320.00 

-166.35 

2.93 


73 

360 . 00 

-194. 12 

4.77 


74 

360. 00 

-194. 12 

6.03 


75 

360. 00 

-182.45 

7.50 


76 

360.00 

-170.79 

7.50 


77 

360 . 00 

-159. 12 

6.03 


73 

360.00 

-159. 12 

4.77 


79 

360. 00 

-170. 79 

3. 30 


30 

360 . 00 

-182.45 

O' . .oO 

CON 1 




ZERO 1 

,2,3; 1, 8 




EC 




1 

1.567 





.7176 





0. 1 




4 

1 . 567 




til 

.7178 




6 

. 1 




7 

. 1 




a 

. 1 




SA 




1 

. 05 





.052 




•Ir« 

- 05 




4 

.05 




5 

.05 




6 

.052 




7 

. 05 
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8 .05 

9 . 05 

10 .0675 
CXQT ELD 

E41 

NEECT=1 


2 10 

11 3 

1 2 

•j* 

NSECT^ 

=2 



18 26 

27 

19 1 


NSECT^ 




34 42 

43 

35 1 


NSECT^ 

=4 



50 53 

59 

51 1 

c» 

NBECT^ 

=5 



1 9 16 8 

1 2 


8 16 

15 7 

1 2 

2 

^4SECT^ 

=6 



17 25 

32 

24 1 

2 

24 32 

31 

23 1 


.'nISECT^ 

=7 



33 41 

48 

40 1 


40 48 

47 

39 1 


NSECT^ 

=8 



49 57 

64 

56 1 

3 

56 64 

63 

55 1 

•3 

NMAT=2 



NSECT= 9 



9 

16 

11 

10 

16 

15 

12 

11 

15 

14 

13 

12 

17 

24 

19 

18 

24 

'T'T 

20 

19 



21 

20 

25 


27 

26 


31 

28 

27 

"T -1 

o ± 

30 

29 

28 


40 

35 

34 

40 

39 

36 

35 

39 

38 

37 

36 

41 

48 

43 

42 

48 

47 

44 

43 

47 

46 

45 

44 

49 

56 

51 

50 


i=rcr 

nr '~i 



wJvJ 

Ux:, 

vJ 1 

55 

54 

53 

52 

57 

64 

59 

53 

64 

63 

60 

59 

63 

62 

61 

60 

65 

/ -1. 

67 

66 

72 

71 

68 

67 

71 

70 

69 

63 

73 

80 

75 

74 

SO 

79 

7 6 

75 
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79 78 77 76 

NSECT=10 

1 9 10 2. 1 9 
6 14 13 5 1 9 

E23 
NMAT=1 
NSECT= 1 
19 12 2 7 
6 14 1 2 2 1 
NSECT=2 
22 30 1 2 2 1 

17 25 1 2 2 7 
NSECT=3 

38 46 1 2 2 1 

33 41 1227 
NSECT=4 

2 10 1 2 4 1 
NSECT=5 

18 26 1 2 4 1 
MSECT=6 

34 42 1 2 4 1 
NSECT=7 

54 62 1 3 2 1 

49 57 1 3 2 7 
NSECT=8 

50 58 1 3 4 1 
CXDT EXIT 

EDI ENCOUNTERED. 

/ 
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repeatable part of SPAR 


LIST,FN=REPT 
CXQT E 

RESET G=386. 

CXQT TOPQ 
CXQT EKS 
CXQT K 
CXQT INV 
CXQT EI6 
RESET INIT=8,NREQ=6,HIST=0 
CXQT DCU 

PRINT 1 VIBR EVAL 
CXQT VPRT 
PRINT VIBR MODE 
CXQT AUS 
SYSVEC;UNIT VEC 
1=1; J=i ,80; 1.0 

DEFINE UN=UNIT VEC 
DEFINE DM=DEM DIAB 
OBJF=XTY(UN,DM) 

DEFINE Vl=viBR MODE 1 1 1,6 
XDIS=UNI0N(V1) 

CXQT GSF 
SOURCE=XDIS AUS 
CXQT DCU 
TOC 1 
CXQT EXIT 
EOI ENCOUNTERED. 

/ 


LIST,FN-SECPROP new section properties 

CXQT TAB file created in GRAD 

BC 

1 1.0923 

2 .324 


4 .3 

5 .3 

lb a O 

y a 3 
8 .3 
SA 

1 .05 

2 .05 

3 .05 

4 .05 

5 . 05 

6 . 05 

7 .05 

8 . 05 

9 . 05 

1 0 . 065 
CXQT EXIT 

EOI ENCOUNTERED. 

/ 
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PHWING ONES IX 

^:INPUT NM=1B, NR=0, NC=0, NK=7, ISPLANE=3, C=69.58, 
NL0AD=72, IDCM=1, IAPLT=0, ISTABLE=1 , I6UST=1 , 

BN=.30,.65, NCCEF=5, NC0L=0,6»1, NKK=7*7, 

IC0NSYS=0, NS=0, ISELECT=1, ICSACT=1 , IBL0CK=2, KVAR=3, 
6L=30000. , ISPECT=2, IPRINT=1 , ISQUT=0 , 

V0=Q876.6,RHD0=0. , DV=1.873E-9, NV=0,NG=1, 

N I NT=5 , RFLD W=0 . , RFCUT= 1.0, UU=B876 . 56 , NOEL= 1 , RH0=5 . 7 1 6E-8 , 
IFLPLT=0, IFRFPLT=1$ 

ALL 1 0 
END 

*ACTINP 

^SELECT NM0DES=6 , NCNEW=0, NQC=1 ,2,3,4,5,6* 

•i:PLTSEL t- 


MODES SEVEN TWELVE 

^ INPUT NM=18, NR=0, NC=0, NK=7, ISPLANE=3, 0=69.58, 
NLDAD=72, IDCM=1, IAPLT=0, ISTABLE=1 , IGUST=1 , 

BN=.30,.65, NCDEF=5, NC0L=0,6*1, NKK=7*7, 

IC0NSYS=0, NS=0, ISELECT=1, ICSACT=1 , IBL0CK=2, KVAR=3, 
BL=30000 . , I SPECT=2 , I PR I NT= 1 , I S0UT=0 , 

V0=8876.6,RH00=O. , DV=1.873E-9, NV=0,NG=1, 

N I NT=5 , RFL0W=0 . , RFCUT= 1.0, UU=8876 . 56 , NVEL= 1 , RH0=5 . 7 1 6E-8 , 
I FLPLT=0 , I FRFPLT= 1 
ALL 1 0 
END 

$ACT I NP 

^SELECT NM0DES=6 , NCNEW=0 , N0C=7 ,8,9,10,11,12^ 

^:PLTSEL $ 


MODES THIRTEEN EIGHTEEN 

INPUT NM=1S, NR=0 , NC=0 , NK=7, ISPLANE=3, 0=69.58, 
NL0AD=72 , I DCM= 1 , I APLT=0 , I ST ABLE= 1 , I GUST= 1 , 

BN=.30,.65, N00EF=5, N00L=0,6*1, NKK=7*7 , 

ICDNSYS=0, NS=0, ISELE0T=1, I0SA0T=1 , IBL00K=2, KVAR=3, 

GL=30000. , ISPE0T=2, IPRINT=1 , IS0UT=0 , 

VO=8a76.6,RHOO=0. , DV=1.873E-9, NV=0,NG=1, 

N I NT=5 , RFL0W=0 . , RFCUT= 1.0, UU=S876 . 56 , N VEL= 1 , RH0=5 . 7 1 6E-8 , 
IFLPLT=0, IFRFPLT=li 
ALL 1 O 
END 

•TACTINP T 

:SSELEOT NMDDES=6 , NCNEW=0 , N00= 13,14,15,16,17,18$ 

■TPLTSEL $ 


MODES NINETEEN TWENTYFOUR 

$ INPUT NM=ia, NR=0, N0=0 , NK=7 , ISPLANE=3, 0=69.53, 

NLD AD=72 , I DOM= 1 , I APLT=0 , I STABLE= 1 , I GUST= 1 , 
BN=.30,.65, NC0EF=5, NC0L=0,6*1, NKK=7*7, 

ICDNSYS=0, NS=0, ISELECT=1 , ICSACT=1, IBL0CK=2, KVAR=3, 
GL=30000 . , I SPECT=2 , I PR I NT= 1 , I S0UT=0 , 

V0=8876. 6,RH00=0. , DV=1.873E-9, NV=0,NG=1, 
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N I NT=5 , RFL0W=0 . , RFCUT= 1.0, UU=8S76 . 56 , N0EL= 1 , RH0=5 . 7 1 6E-8 , 

I Fl_PLT=0 , I FRFPLT= 1 
ALL 2 0 
END 

S^ACTINP ^ 

^SELECT NM0DES=6,NCNEW=0,N0C=1 ,2,3,4,5,6$: 

:^PLTSEL ^ 

O 

-2 

MODES TWENTYFIVE THIRTY 

INPUT NM=18, NR=0, NC=0, NK=7, ISPLANE=3, 0=69.58, 

NLDAD=72 , I DCM= 1 , I APLT=0 , I STABLE= 1 , I GUST= 1 , 

BN=.30,.65, NC0EF=5, NC0L=0,6*1, NKK=7*7, 

IC0NSYS=0, NS=0, ISELECT=1, ICSACT=1 , IBL0CK=2, KVAR=3, 

8L=30000 . , I SPECT=2 , I PR I NT= 1 , I S0UT=0 , 

V0=8876.6,RHQ0=0. , DV=1.873E-9, NV=0,NG=1, 

N I NT=5 , RFL0W=0 . , RFCUT= 1.0, UU=8876 . 56 , NVEL= 1 , RH0=5 . 7 1 6E-B , 

I FLPLT=0 , I FRFPLT= 1 
ALL 2 0 
END 

:#ACTINP ■$ 

^:SELECT NM0DES=6 , NCNEW=0 , NDC=7 ,8,9,10,11,12-4: 

$PLTSEL 

O 

—2 

MODES THIRTYONE THIRTYSIX 

INPUT NM=18, NR=0, NC=0, NK=7, ISPLANE=3, 0=69.58, 

NL0AD=72, IDCM=1, IAPLT=0, ISTABLE=1 , IGUST=1 , 

BN=.30,.65, N00EF=5, N00L=0,6*1, NKK=7*7, 

I00NSYS=0, NS=0, ISELECT=1, ICSACT=1, IBL00K=2, KVAR=3, 

GL=30000. , ISPECT=2, IPRINT=1 , IS0UT=O , 

V0=B876 . 6 , RHO0=O . , D V= 1 , B73E-9 , N V=0 , NG= 1 , 

N I NT=5 , RFL0W=0 . , RFCUT= 1.0, UU=8876 . 56 , NVEL= 1 , RH0=5 . 7 1 6E-8 , 

IFLPLT=0, IFRFPLT=14: 

ALL 2 O 
END 

$ACTINP t- 

4:SELECT NM0DES=6 , NCNEW=0 , NOC= 13,14,15,16,17,184: 

^PLTSEL 4: 

0 

-2 

MODES THIRTYSEVEN FDRTYTWO 
4: INPUT NM=ia, NR=0, NC=0, NK=7 , ISPLANE=3, 0=69.58, 

NLO AD=72 , I DCM= 1 , I APLT=0 , I STABLE= 1 , I GUST= 1 , 

BN=.30,.65, NC0EF=5, NC0L=0,6*1, NKK=7*7 , 

IC0NSY3=0, NS=0, ISELECT=1 , ICSACT=1, IBL0CK=2, KVAR=3, 

GL=30000. , ISPECT=2, IPRINT=1 , IS0UT=0 , 

V 0 =B 876 . 6 ,RHa 0 = 0 . , DV=l.S73E-9, NV=0,NG=1, 

N I NT=5 , RFLQW=0 . , RFCUT= 1.0, UU=S376 . 56 , NVEL= 1 , RH0=5 . 7 1 6E-B , 

IFLPLT=0, IFRFPLT=14: 

ALL 30 * 

END 

4:ACTINP 4: 

4SELECT NM0DES=6 , NCNEW=0 , N0C=1 ,2,3,4,5,64: 

4:PLTSEL 4: 
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MODES FORTYTHREE FORTYEIGHT 

INPUT NM=ia, NR=0, NC=0, NK=7, ISPLANE=3, C=69-58, 

NLQ AD=72 , I DCM= 1 , I APLT=0 , I ST ABLE= 1 , I 6UST= 1 , 

BN=.30,.65, NC0EF=5, NC0L=0,6*1, NKK=7*7, 

IC0NSYS=0, NS=0, ISELECT^l, ICSACT=1, IBL0CK=2, KVAR=3, 
GL=30000 . , I SPECT=2 , I PR I NT= 1 , I S0UT=0 , 

V0=9876.6,RH00=0. , DV=1.873E-9, NV=0,NG=1, 

N I NT=5 , RFL0W=0 . , RFCUT= 1.0, 88=8876 - 56 , NVEL= 1 , RHD=5 . 7 1 6E~S , 
IFLPLT=0, IFRFPLT=1# 

ALL 3 O 
END 

i ACT IMP ^ 

^SELECT NM0DES=6 , NCNEW=0 , N0C=7 ,8,9,10,11,12$: 
tPLTSEL $• 


MODES FORTYNINE FIFTYFOUR 

INPUT NM=18, NR=0, NC=0, NK=7, ISPLANE=3, C=69.58, 
NL0AD=72, IDCM=1, IAPLT=0, ISTABLE=1 , IGUST=1 , 

BN=.30,.65, NC0EF=5, NC0L=0,6*1, NKK=7*7, 

IC0NSYS=0, NS=0, ISELECT=1, ICSACT=1, IBL0CK=2, KVAR=3, 
GL=30000 . , I SPECT=2 , I PR I NT= 1 , I S0UT=0 , 

V0=8876 . 6 , RH00=O . , DV= 1 . 873E-9 , NV=0 , NG= 1 , 

N I NT=5 , RFLO W=0 . , RFCUT= 1.0,88=8876.56, N VEL= 1 , RH0=5 . 7 1 6E-8 , 
IFLPLT=0, IFRFPLT=1:$ 

ALL 3 0 
END 

$ACTINP ^ 

^SELECT NM0DES=6 , NCNEW=0 , NOC= 13,14,15,16,17,18^ 

$PLTSEL 

0 

12 

EOI ENCOUNTERED. 
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LIST,FN=TAFE31 
PHWING TEST MODEL 
ANYTH I NB 
<3E15.8) 

$DL I NPT NSECTNS= 1 , NSECTNA= 1 , NDF=3 , NN0DES=40 , NN0DET=40 , 
IDLAT=1 ,NELAST=1S, IPLTMS=0, 

IPRINT=0, ISYM=+1 , 

NPMX=40, IPLTB0X=0, IPLTN0D=0 ,SCX=30. ,SCY=30. , IDCM=1 , ITOC=l^ 
DISP -1 1 
MODES 1 1 
HHD 1 1 
SPLINE 1 1 
END 

1 40 1 0. 0. 1. 

1 48 1 


67-5000 

0.0000 

1 

o 

-5 

4 

37.5000 

0.0000 

2 

3 

-5 

4 

7.5000 

0. 0000 

3 


*~5 

4 

-22. 5000 

0.0000 

4 


“5 

4 

81.5700 

40.0000 

5 

o 

-5 

4 

53. 6100 

40. 0000 

6 

3 

-5 

4 

25.6400 

40.0000 

7 

3 

-5 

4 

-2.3200 

40.0000 

8 

3 

—5 

4 

95.6400 

80.0000 

9 

3 

““5 

4 

69.7100 

80.0000 

10 

3 

--5 

4 

43.7900 

80.0000 

11 

3 

--5 

4 

17.8600 

80.0000 

12 

3 

—5 

4 

109.7100 

120.0000 

13 

3 

-5 

4 

85.8200 

120.0000 

14 

-5** 

o 

--5 

4 

61.9300 

120.0000 

15 

3 

-5 

4 

38.0400 

120.0000 

16 

3 

•“5 

4 

123.7800 

160. 0000 

17 

T 

O 

“•5 

4 

101.9200 

160.0000 

18 


-5 

4 

80.0700 

160. 0000 

19 


—5 

4 

53.2200 

1 60 . 0000 

20 


—5 

4 

137. 8400 

200. OOOO 

21 

3 

*-5 

4 

118.0300 

200.0000 

22 

-r 

-5 

4 

98.2200 

200 . 0000 

^*T»' 

*T 

~-5 

4 

78. 4000 

200.0000 

24 


-5 

4 

151 . 9100 

240 . 0000 

25 


~5 

4 

134. 1400 

240. 0000 

26 

3 

-5 

4 

1 16.3600 

240. OOOO 

27 

T 

—5 

4 

98 . 5800 

240. OOOO 

28 


-5 

4 

165. 9800 

280. OOOO 

29 

3 

*“5 

4 

150. 2400 

230. OOOO 

30 


•**5 

4 

1 34 . 5000 

280. OOOO 

31 

y. 

~5 

4 

113. 7600 

230. OOOO 


o 

^5 

4 

180- 0500 

320. OOOO 



-“5 

4 

166. 3500 

320 . OOOO 

34 


~5 

4 

152. 6400 

320 . OOOO 

35 

•— » 

cr 

4 

138. 9400 

320.0000 

36 

♦J> 

-^5 

4 

194. 1200 

360 . OOOO 

37 

3 

“~5 

4 

182.4500 

360.0000 

38 

3 

~5 

4 

170.7900 

360. OOOO 

39 

3 

—5 

4 

159- 1200 

360. OOOO 

40 

•u> 

—5 

4 


EG I ENCOUNTERED. 

/ 


A- 30 



LIST,FN=TAPE32 
PHWING TEST MODEL 
ANYTHING 
(3E15.G) 

I NPT NSECTNS= 1 , NSECTNA= 1 , NDF=3 , NNQDES=40 , NN0DET=40 , 

I DL AT= 1 , NELAST=6 , I PLTMS=0 , 

IPRINT=0, ISYM=+1 , 

NPMX=40, IPLTB0X=0, IPLTN0D=0,SCX=30. ,SCY=30. , IDCM==1 , ITOC=l:$ 
DISP -12 
MODES 2 2 
HHD 2 2 
SPLINE 2 2 
END 


67.5000 

0 . 0000 

1 

3 

-5 

4 

37.5000 

0. 0000 

2 

3 

—5 

4 

7. 5000 

0. 0000 

3 


“5 

4 

-22. 5000 

0.0000 

4 


-5 

4 

81.5700 

40.0000 

5 

“T 

-”5 

4 

53.6100 

40.0000 

6 

3 

~5 

4 

25.6400 

40.0000 

7 


“5 

4 

-2.3200 

40.0000 

8 

3 

—5 

4 

95. 6400 

80 . 0000 

9 

3 

“5 

4 

69.7100 

80. 0000 

10 


“*5 

4 

43. 7900 

80. 0000 

11 


“5 

4 

17.3600 

SO. 0000 

12 


““5 

4 

109.7100 

120.0000 

13 

o 

“5 

4 

85. 8200 

120.0000 

14 


~5 

4 

61.9300 

1 20 . 0000 

15 

T 

-5 

4 

38. 0400 

120. 0000 

16 

-T 

■_« 

— 5 

4 

123.7800 

160. 0000 

17 

3 


4 

101. 9200 

160. 0000 

18 


—5 

4 

80 . 0700 

1 60 . 0000 

19 

3 

— 5 

4 

58.2200 

160.0000 

20 

o 

“5 

4 

1 37 . 8400 

200. 0000 

21 

o 

—5 

4 

118.0300 

200.0000 



-5 

4 

98.2200 

200.0000 


T 

“•5 

4 

78. 4000 

200 . 0000 

24 

3 

“■5 

4 

151 .9100 

240 . 0000 

25 


“5 

4 

134. 1400 

240. 0000 

26 

3 

“5 

4 

116. 3600 

240.0000 

27 

3 

“5 

4 

98.5800 

240. 0000 

23 


”“5 

4 

165.9800 

280.0000 

oa 

X. / 

-j* 

“5 

4 

150.2400 

280.0000 

30 


—5 

4 

134. 5000 

280. OOOO 

31 

3 

“5 

4 

118. 7600 

280. 0000 


3 


4 

180. 0500 

320. 0000 

"TT 

3 

“"5 

4 

1 66 . 3500 

320. 0000 

34 

3 

—5 

4 

152,6400 

320. 0000 

35 

3 

“5 

4 

138.9400 

320. 0000 

36 

3 

“5 

4 

194. 1200 

360. 0000 

37 


~5 

4 

182. 4500 

360 . 0000 

3B 


”*5 

4 

170. 7900 

360. 0000 

39 

3 

—5 

4 

159.1200 360.0000 

ENCOUNTERED. 

40 

o 

“"5 

4 
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LIST,FN=TAPE33 
PHWING TEST MODEL 
ANYTHING 
(3E15.8) 

«=DL I NPT NSECTNS= 1 , NSECTNA= 1 , NDF=3 , NN0DES=40 , NN0DET=40 
IDLAT=1 ,NELAST=6, IPLTMS=0, 

IPRINT=0, ISYM=+1 , 

NPMX=40,IPLTB0X=0,IPLTN0D=0,SCX=30. ,SCY=30. , IDCM=1 , ITOC=l^ 
DISP -1 3 
MODES 3 3 
HHD 3 3 
SPLINE 3 3 
END 

I 40 1 0. O. 1. 

1 48 1 


67.5000 

0.0000 

1 

-r 

O 

-5 

4 

37.5000 

0.0000 

2 

w 

-5 

4 

7.5000 

0.0000 

o 

3 

*~5 

4 

-22.5000 

0. 0000 

4 

-?r 

O 

—5 

4 

81.5700 

40.0000 

5 


—5 

4 

53.6100 

40.0000 

6 

“T 

-5 

4 

25.6400 

40. 0000 

7 

3 

^5 

4 

-2.3200 

40.0000 

8 


-5 

4 

95.6400 

80. 0000 

9 

'T 

-5 

4 

69.7100 

80. 0000 

10 


—5 

4 

43.7900 

80.0000 

1 1 


—5 

4 

17.8600 

80.0000 

12 

3 

-5 

4 

109.7100 

120.0000 

13 

o 

*~5 

4 

85.8200 

120.0000 

14 

•u> 

--5 

4 

61.9300 

120.0000 

15 

3 

~*5 

4 

38.0400 

120. 0000 

16 

c* 

““5 

4 

123.7800 

1 60 . OOOO 

17 

3 

~5 

4 

101.9200 

160. 0000 

18 


“*5 

4 

80 . 0700 

160. 0000 

19 


“**5 

4 

58.2200 

160. 0000 

20 

■u» 

"“5 

4 

137. 8400 

200. 0000 

21 

3 

“"5 

4 

118. 0300 

200. 0000 

22 


“•5 

4 

98. 2200 

200. 0000 



~5 

4 

78. 4000 

200. 0000 

24 


~5 

4 

151 .9100 

240. 0000 

25 


-5 

4 

134. 1400 

240. 0000 

26 

o 

-5 

4 

116. 3600 

240. 0000 

27 


““5 

4 

93. 5800 

240 . 0000 

23 

3 

~5 

4 

165. 9300 

280.0000 

29 

3 

~5 

4 

150. 2400 

280 . OOOO 

30 

-T 

—5 

4 

134.5000 

280 . OOOO 

31 

3 

~5 

4 

1 18.7600 

280. OOOO 

"T'T* 

•j> 

~*5 

4 

1 SO . 0500 

320. OOOO 

33 


-5 

4 

166. 3500 

320. OOOO 

34 

3 

*"5 

4 

152.6400 

320.0000 

"Tcr 

7:; 

cr 

~vJ 

4 

133. 9400 

320. OOOO 

36 

•.!> 

“■5 

4 

194. 1200 

360. OOOO 

37 


5 

4 

182.4500 

360. OOOO 

33 

3 

—5 

4 

170.7900 

360. OOOO 

39 

3 


4 

159.1 200 

360. OOOO 

40 

0 

—5 

4 


EDI ENCOUNTERED. 
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LIST,FN=TAPE5T 

PHWING TEST MODEL _ , 

AT I NP ACAP=22500 . O , REFCHD=69 . 44 , 

FMACH=.700, 

NDELT=1 ,NP=1 ,NB=0,NRF=7,RFREQ=0. 0,0. 05,0. 15,0.30,0.5,0.75,1.05, 
NSTR I P= 1 2 , NS V=0 , NBV=0 , NMODES= 1 8 , ND AT A=0 , 

NELAST= 1 8 , REFSFN= 1.0, NB0X=48 , NSNC=65 , 

X ZER0=-22 . 5 , NR=0 , 

NC=0 , 1 DCM= 1 , I TOC-- 1 
HHD 1 O 
AERO 1 1 
AIC 0 1 
END 

-22.5 67.5 159.12 194.12 0.0 360.0 

7.02 69.51 13 5 1.0 
0.0 0.25 0.5 0.75 1.0 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.666 0.732 0.798 0.864 0.93 1.0 
1 4 5 8 9 12 13 16 17 20 21 24 
25 28 29 32 33 36 37 40 41 44 45 48 
SECOND CASE 

#DL AT I NP ACAP=22500 . 0 , REFCHD=69 . 44 , 

FMACH=.700, 

NDELT= 1 , NP= 1 , NB=0 , NRF=7 , RFREQ=0 .0,0.05,0.15,0.30,0.5,0.75,1.05, 
NSTR I P= 1 2 , NSV=0 , NBV=0 , NMODES= 1 8 , NDATA=0 , 

NEL AST= 1 8 , REFSPN= 1.0, NB0X=4S , NSNC=65 , 

X ZER0=-22. 5 , NR=0 , 

NC=0 , 1 DCM= 1 , 1 T0C=2 t- 
HHD 2 0 
AERO 2 2 
AIC 1 O 
END 

-22.5 67.5 159.12 194.12 0.0 360.0 

7.02 69.51 13 5 1.0 
0.0 0.25 0.5 0.75 1.0 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.666 0.732 0.798 0.864 0.93 1.0 
1 4 5 a 9 12 13 16 17 20 21 24 
25 28 29 32 33 36 37 40 41 44 45 48 
THIRD CASE 

i^DLATINP ACAP=22500.0,REFCHD=69.44, 

FMACH=.700, 

NDELT= 1 , NP= 1 , NB=0 , NRF=7 , RFREQ=0 .0,0.05,0.15,0.30,0.5,0.75,1.05, 
NSTR I P= 1 2 , NSV=0 , NBV=0 , NMODES= 1 3 , NDATA=0 , 

NELAST= 1 8 ‘ REFSPN= 1.0, NB0X=--48 , NBNC=65 , 

XZER0=-22.5,NR=0, 

NC=0,IDCM=1 ,IT0C=3 f 
HHD 3 O 
AERO 3 3 
AIC 1 O 
END 

-22.5 67.5 159.12 194.12 0.0 360.0 

7.02 69.51 13 5 1.0 
0.0 0.25 0.5 0.75 1.0 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.666 0.732 0.798 0.864 0.93 1.0 
1 4 5 8 9 12 13 16 17 20 21 24 
25 28 29 32 33 36 37 40 41 44 45 48 
EDI ENCOUNTERED. 

/ 
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